Electric-Field-Induced Lifting of the Valley Degeneracy in 
a-(BEDT-TTF) 2 I 3 Dirac-Like Landau Levels 
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The relativistic Landau levels in the layered organic material a— (BEDT-TTF)2l3 [BEDT- 
TTF=bis(ethylenedithio)tetrathiafulvalene] are sensitive to the tilt of the Dirac cones, which, as 
in the case of graphene, determine the low-energy electronic properties under appropriate pressure. 
We show that an applied inplane electric field, which happens to be in competition with the tilt of 
the cones, lifts the twofold valley degeneracy due to a different level spacing. The scenario may be 
tested in infrared transmission spectroscopy. 
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The study of relativistic effects in condensed matter 
systems has become a major issue since the discovery 
of a particular quantum Hall effect in graphene [j], Q, 
which may be understood in terms of relativistic two- 
dimensional (2D) charge carriers at the Fermi level 
These relativistic carriers occur in the vicinity of the crys- 
tallographic points K and K' within the first Brillouin 
zone (BZ), where the valence band touches the conduc- 
tion band. The electronic energy dispersion is then linear 
in the wave vector k, e(fc) = hvpk, in terms of the Fermi 
velocity vf ~ 10 6 m/s, which is roughly 300 times less 
than the velocity of light. 

Another system, less known than graphene, where 
2D relativistic charge carriers are expected to de- 
scribe the low-energy electronic properties, is the or- 
ganic layered compound a— (BEDT-TTF) 2 l3 [BEDT- 
TTF=bis(ethylenedithio)tetrathiafulvalene] 0, H, 0. 
The electronic structure is, however, more involved than 
in graphene due to the presence of four sites per unit 
cell with highly anisotropic overlap integrals. An addi- 
tional complication stems from the presence of strong 
electronic correlations in these materials, which lead to 
insulating phases with modulated charge density at low 
temperatures [7|. Fermi-liquid- type physics may be in- 
duced when pressure is applied to this 3/4-filled system, 
in which case the Fermi level resides at the contact points 
of the upper two bands, which form Dirac cones 0, [||. 
Correlation effects, as well as different on-site energies, 
may then be taken into account within a mean-field the- 
ory and renormalize the band parameters 0]. 

A particularly interesting feature of Dirac cones in 
a— (BEDT-TTF)2l3, which occur in two distinct points 
within the first BZ, is their tilt in the momentum-energy 
space. This tilt may be understood as due to an inter- 
play between a highly anisotropic nearest-neighbor hop- 
ping, which drives the Dirac points away from points 
of high crystallographic symmetry, and a next-nearest- 
neighbor hopping which is of comparable size as the 
nearest- neighbor one 9]. In a magnetic field, the spec- 
trum is organized in relativistic Landau levels (LL) as 
in graphene, but the tilt of the Dirac cones reduces the 



effective Fermi velocity and, thus, the LL spacing [9I. fTol] . 

In spite of this theoretical understanding and in con- 
trast to graphene, strong experimental evidence for 
(tilted) Dirac cones in a-(BEDT-TTF) 2 I 3 under pres- 
sure is yet lacking. An unusual T 2 dependence of the 
carrier density, observed in transport measurements ll| , 
may be interpreted as an indication of relativistic car- 
riers 0, H|. This temperature dependence is in com- 
petition with a T~ 2 dependence of the mobility, such 
that the electric conductivity is roughly temperature- 
independent. Very recently, magneto-transport measure- 
ments have revealed the existence of a zero-energy LL in 
a-(BEDT-TTF) 2 I 3 [H, in agreement with theoretical 
predictions [9(. 

Another means of reducing the LL spacing (in 
graphene), different from the above-mentioned tilt of 
the Dirac cones due to a particular combination of 
the hopping parameters, has been proposed by Lukose, 
Baskaran, and Shankar [l3j |. The application of an in- 
plane electric field reduces indeed the LL spacing by a fac- 
tor [1 - (£/v F B) 2 ] 3 / 4 . This effect is a direct consequence 
of the covariance of the Dirac equation for fermions in an 
electromagnetic field and may be understood in terms of 
a Lorentz transformation into a reference frame, where 
the electric field vanishes. 

Here, we show that the combination of an inplane elec- 
tric field and a perpendicular magnetic field B = Be z 
may be viewed as an effective change of the tilt of the 
Dirac cones. Indeed, it renormalizcs the tilt parameter, 
which becomes not only dependent on the electric field 
but also on the valley index £ = ±1. The renormalized 
tilt parameter affects the effective cyclotron energy and 
the LL spectrum. The latter depends explicitly on the 
valley index, and, thus, the inplane electric field lifts the 
twofold valley degeneracy for the LLs n 7^ 0. The cen- 
tral LL n — remains at zero energy in both valleys 
and maintains its twofold valley degeneracy. A quantum 
Hall measurement, which one may first think of to reveal 
the lifted valley degeneracy, is a difficult task due to the 
lack of single-layer a-(BEDT-TTF) 2 I 3 . An alternative 
and more realistic experimental test of the proposed see- 
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nario may be infrared transmission spectroscopy, which 
has been proven to be a powerful tool in the determina- 
tion of relativistic LLs in graphene 14, 

In the absence of a magnetic field, the low-energy elec- 
tronic properties of a— (BEDT-TTF)2l3 with weak layer 
coupling may be described by the minimal form of the 
generalized 2D Weyl Hamiltonian for the valley £ 



Hf=i (w ■ qcr° 



w x q x a 



w, 



(1) 



in terms of the 2x2 Pauli matrices a x , a v , and cr° = 1 
(we use a system of units with U = 1). The velocities 
wo = (wo x ,WQy) and w = (w x ,w y ) may be viewed as 
effective parameters. Estimates, calculated within the 
model proposed by Hotta and based on the overlap 
integrals of Ref. |l7j], yield w x ~ 2.14 eVA, Wy ~ 0.22 
eVA, w Qx ~ -0.0075 eVA, and w 0y ^ 0.074 eVA @. 

In order to account for the perpendicular magnetic field 
Be z = V x A, one may use the Peierls substitution 



n = 



eA, 



(2) 



which is a valid approximation when the lattice spacing 
[~ 1 nm in a— (BEDT-TTF)2l3] is much smaller than the 
magnetic length l B = 1/s/eB ~ 26/VB[T] nm. With 
the help of the ladder operators 



: (W X U X - iWylly) , 



l x Wy 



a f = 7= (w x IL Q 
yJ2w x w y 

one obtains the Hamiltonian Q 
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where 
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in terms of the effective tilt parameter 
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The angle ip denotes the angle between the x-axis and the 
direction of the effective tilt {wo x /w x ,w 0y /w y ), renor- 
malized by the Fermi velocities w x and Wy. in the x- and 
y-direction, respectively. A semiclassical [9( as well as a 
full quantum treatment [l(| of the Hamiltonian ((4| yields 
the energy spectrum 



in = A V^ ( l_^)3/4 



2n. 



(6) 



where A = ± indicates the band index. The spectrum, 
which has the same structure as that of relativistic elec- 
trons in graphene, does not depend on the valley index 



£ and is, thus, twofold (fourfold, when one takes into ac- 
count the electron spin) degenerate, in addition to the 
macroscopic LL degeneracy characterized by the density 
of flux quanta tib = 1/2ttI b . 

In the presence of an inplane electric field E = Ee^ in 
the direction e|| = (cos 9, sin 6) in the xy-plane, one needs 
to add a term eE ■ ra° — eEx 1 cr° to the Hamiltonian 
(j4|), where x' = xcosd + ysm9. For non-tilted Dirac 
cones (wq = 0) and for an isotropic Fermi velocity vp — 
w x = w y , as in the case of graphene, the electric field 
leads to the same energy spectrum ([B]) . One simply needs 
to replace the tilt parameter wq — > E/vfB [131 ] and to 
add a term (E/B)k, where k is the wave vector in the 
direction = (— sin#, cos 9) perpendicular to e||. Here, 
one naturally chooses the Landau gauge A = Bx'e±, 
such that the full Hamiltonian remains translationally 
invariant in the y'-direction (ej_) and the associated wave 
vector k is, therefore, a good quantum number. 

An inplane electric field in a-(BEDT-TTF) 2 I 3 , i.e. 
when the Fermi velocity is anisotropic and the Dirac 
cones are tilted, is slightly more involved. The compli- 
cation stems from the competition of three characteristic 
directions: the main axes of the anisotropic Fermi veloc- 
ities w x and w y , which define the reference frame; the 
direction of the tilt described by <p\ and the direction of 
the applied electric field given by the angle 9. One may 



then express the variable x' — kl B 



IL>,a+) 



-\Jwyjw x (a + aT 



H y (a, a))l B , where 
9 + iyjw x /w y (a - 



aJ) cos9] is the y'-component of the gauge-invariant mo- 
mentum ([2]) . The final Hamiltonian in the presence of an 
electric field reads 

E. 



H, 



E/B 



B 



(7) 



where H'^ is the same as that of Eq. (0| if one replaces 
the tilt parameter wq exjp(iip) by 



Wx 



Here, the renormalized tilt velocity is given by 

ExB 



B 2 



(8) 



(9) 



and the angle (p^ is related to the angle between the renor- 
malized tilt velocity and the a;-axis. 

Eq. (J9j) is the central result of this paper. The ad- 
ditional term in Eq. ([9]) is nothing other than the drift 
velocity V£> of a charged particle in the presence of both 
a magnetic and an electric field This drift velocity 
may therefore, on its own, tilt the original Dirac cones 
in the same manner as the velocity Wo, which is due to 
the special combination of the nearest-neighbor and next- 
nearest-neighbor overlap integrals in the tight-binding 
model. The intervening electric field depends effectively 
on the valley index £, and so does the LL spectrum, 

* s (E) = X^[l-ME) 2 f A ^ ^ 



Ib 



<2n 



B 



k. (10) 
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One may distinguish three different instances for the 
LL spectrum. 

(i) In the case with no original tilt (w = 0), such 
as for unconstrained graphene, one recovers the valley- 
independent result of Ref . [l3| , where the decrease of the 
effective Fermi velocity is due to the Lorentz invariance of 
the system. Indeed, the effective magnetic field is reduced 
by a factor 



y/1 - (E/v F B) 2 



(11) 



as a consequence of a Lorentz transformation into the 
reference frame, which moves at the drift velocity and 
where the electric field vanishes. Because the energy is 
proportional to 1/Ib oc \f~B~, this leads to a reduction 
by a factor [1 — (E/vfB) 2 } 1 / 4 of the LL spacing. How- 
ever, the spectrum needs to be multiplied by another 
factor (fTT|) when measured in the lab frame, due to the 
Lorentz transformation of the energy. The total reduc- 
tion is, therefore, given by the factor [1 — (E/vfB) 2 ] 3 / 4 , 
in agreement with Eqs. © and (fTTH) . The LL spectrum 
collapses once the drift velocity vd — E/B exceeds the 
Fermi velocity Vp, which is an upper bound in analogy 
with the velocity of light. Physically, the system enters 
then into a regime which is dominated by the electric 
field, and where a Lorentz transformation ma y b e found 
into a reference frame with no magnetic field [18[ . 

(ii) In the case of originally tilted Dirac cones (wo 7^ 
0) and no electric field, one recovers the LL spectrum 
discussed in Refs. [§] and [l(J. As for the above case (i), 
the spectrum is valley-degenerate and collapses beyond 
a critical tilt value given in terms of the tilt parameter 
Wq = 1. Indeed, for w < 1 the semiclassical trajectories 
are closed orbits (ellipses), whereas those for w a > 1 are 
open parabolas, which would give rise to a continuous 
spectrum. The two regimes of closed and open orbits 
may therefore be identified with the relativistic ones of 
the case (i) dominated by the magnetic and the electric 
field, respectively. 

(Hi) The most interesting situation arises when both 
the electric field and tilted Dirac cones are present, such 
as in a-(BEDT-TTF) 2 I 3 . In this case, the valley de- 
generacy is lifted, and the energy spectra are different 
for the two valleys. The effect is most pronounced when 
w _L E, in which case the electric field redresses the cone 
in one valley (say £ = + ^23| ) while enhancing the tilt in 
the other one, £ = — , (Fig. [1]). This leads to a larger LL 
spacing in £ = +, whereas that in £ = — decreases, and 
the critical electric field, beyond which the LL spectrum 
collapses and becomes continuous, is different in the two 
different valleys. Notice that the zero-energy LL remains 
unaffected by the electric field and maintains, therefore, 
its original twofold valley degeneracy. 

A particularly interesting value of the electric field is 
E C omp = l w o|^j when the drift velocity exactly com- 
pensates wo. Notice that because wq is almost paral- 
lel to the y-axis (ip ~ 91°), the electric field is then to 
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FIG. 1: Sketch of the valley-dependent LL spectrum for 
tilted Dirac cones in the presence of an electric field (thick 
lines). We omitted the inclination of the LLs due to the 
term (E/B)k, which lifts the LL degeneracy. The dashed 
lines schematically represent the tilted cones in the two val- 
leys £ = + and — in the absence of an electric field. The 
cones in the two valleys are tilted in opposite directions in 
the momentum-energy space, whereas the electric field acts 
in the same direction. The LL spectrum in the presence of 
an electric field in £ = + (£ = — ) is that of a cone with a 
decreased (increased) tilt (full lines). 



be chosen in the x-direction. In this case, the LL spac- 
ing is maximal in one valley, with a cyclotron energy of 



yj2w x w y /Ib, whereas it is e r 



y/2W X Wy[l 



4u>§] 3 / 4 /Ib in the other one. With the above-mentioned 
values of the velocities, one has |wo| ~ 10 5 m/s, and the 
compensation condition is, therefore, fulfilled for an elec- 
tric field of E comp ~ 10 6 V/m at a typical magnetic field 
of 10 T. Inplane electric fields of this strength (and larger) 
have indeed been obtained in graphene [191 ] , and one may 
expect that they are in the accessible experimental range 
also in a— (BEDT-TTF)2l3, which is expected to have 
similar screening properties as graphene. A rough esti- 
mate of the tilt parameter in a— (BEDT-TTF) 2 l3 yields 
w Q ~ 0.33 |9j, and one may therefore expect a rather 
large ratio e m in/^max ~ 0.6, i.e. a 40% modulation. 

As an experimental test of the proposed scenario, one 
may first think of a quantum Hall measurement, where 
one would expect that the plateau series at the filling fac- 
tors v — n e i/riB = ±2(2n+l) of the relativistic quantum 
Hall effect evolves into a series of v = ±2(n + 1) once a 
sufficiently large electric field is applied. Here, n e i is the 
2D carrier density, and we have taken into account the 
twofold spin degeneracy. Such a measurement would, 
however, require a single layer of a— (BEDT-TTF)2l3, 
which has not been fabricated yet. Another more realistic 
experimental test of the proposed scenario consists of an 
infrared spectroscopic measurement, which has been per- 
formed with great succes in epitaxial 14j and exfoliated 
graphene 15| and which reveals the relativistic LL struc- 
ture of these materials. The dipolar coupling is sensitive 
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to the transitions A, n — > A', nil, and one may therefore 
access both intraband (A' = A) and interband (A' = —A) 
transitions. In the case of 3/4-filled a-(BEDT-TTF) 2 I 3 , 
where the n = LL is half-filled, one would expect ab- 
sorption lines at energies (for all integers n) 

Al(E) = [l - ^{Eff* (V^+V^Tl), 

(12) 

which split into two lines each when an inplane electric 
field is applied. A linear expansion in E of Eq. (|12p yields 
a relative splitting r/ = [A+(E) - A^(E)]/A^(E = 0) ~ 
E/vpB of the lines if the direction of E is appropriately 
chosen (close to the e x direction, where the tilt is maxi- 
mal) and if the LL broadening is sufficiently small. Here, 
v F = y/w x w y [l ~ w^(E) 2 } 3 / 4 ~ 10 5 m/s is the renormal- 
ized Fermi velocity which accounts for the tilt as well 
as for the anisotropy of the Dirac cones in a-(BEDT- 
TTF) 2 I 3 . 

Notice, however, that the strength of the inplane elec- 
tric field is delimited by breakdown processes, which 
may rearrange the LL filling [2(j. Above a critical elec- 
tric field E c , once the associated potential energy ac- 
quired over the characteristic length ^ is a substan- 
tial fraction 7 of the inter-LL spacing 5 n = e n +\ — e n , 
eE c ls = 7<5n, an electron may tunnel from one LL n 
to a neighboring one n + 1. In GaAs heterostructures, 
7 ~ 0.1 is a reasonable estimate to account for the ob- 
served breakdown [21| . If one considers a similar limi- 
tation in a— (BEDT-TTF)2l3, the breakdown occurs at 
E c ~ jvpB ~ 10 5 V/m < Ecomp, and the maximal 
relative strength of the valley-degeneracy lifting may be 
estimated from (fT2j) as 7, which yields a ~ 10% effect, 
for a cyclotron energy of roughly e ~ ZAyJ B\T] meV (10 
meV at S = 10 T). 

We finally point out that a similar valley-degeneracy 
lifting may occur un graphene, where a tilt of the Dirac 
cones may, in principle, be induced under uniaxial strain 
0. The effective tilt parameter has beeen estimated as 
wo ~ 0.6e, where e is the relative variation of the bond 
length under uniaxial strain [9], which may be on the 
order of 10-20% before the graphene sheet cracks 
The order of magnitude for lifting the valley degeneracy 
is -q ~ (E / Bvp)e, and if one assumes a breakdown pa- 
rameter 7 ~ 0.1, this yields a ~ 1% effect for a critical 
electric field of roughly 10 6 V/m, which is roughly ten 
times larger than that of a-(BEDT-TTF) 2 I 3 , due to the 
larger Fermi velocity in graphene. 

In conclusion, we have shown that a simultaneous tilt 
of the Dirac cones and an applied electric field in the 
2D plane is capable of lifting the twofold valley degen- 
eracy of the LL spectrum. Whereas one may induce a 
tilt of the Dirac cones in graphene by application of a 
uniaxial strain those in a— (BEDT-TTF) 2 I 3 are nat- 
urally tilted under moderate pressure, due to a particular 
interplay between the nearest-neighbor and next-nearest- 



neighbor hopping integrals 0, H, @, Q . The inplane elec- 
tric field acts, in the presence of a perpendicular mag- 
netic field, as an additional effective tilt, which is valley- 
dependent and which competes with the original tilt. If 
the direction of the inplane field is properly chosen, it 
redresses the Dirac cone in one valley £ while enhancing 
the tilt in the other one — £. The LL spectrum is there- 
fore stretched in £ and compressed in — £. The effect may 
be observable in infrared transmission spectroscopy. 
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